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ABSTRACT 

This work was undertaken to determine the pressure drops in the 
oxidizer flow passages of a compact dual -combust ion liquid rocket 
engine by means of fluid dynamic model tests, using air instead of 
liquids, to establish design data. Investigated also were flow dis- 
tributions and flow instabilities by measurements and flow visuali- 
zation methods. Changes are indicated to improve the original design. 

The air tests were carried out over a range of Reynolds numbers 
to establish the effect of compressibility on the pressure drop loss 
coefficients. Water tests were undertaken to check the validity of 
the results and of the proposed calculating method. 

The tests were conducted at the Turbo-Propulsion Laboratory of 
the Naval Postgraduate School, Monterey, California. 
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SECTION 1 



INTRODUCTION 

The science of astronautics has progressed at an exponential rate 
in the last decade. The development of liquid rocket engines capable 
of producing large thrusts at high specific impulse has been responsi- 
ble for this progress. These engines require very large flow rates of 
fuel and oxidizer. To keep their size and weight at a minimum the com- 
bustion pressures are high, limited only by structural and turbo-pump 
design considerations. The flow passages inside the units are intri- 
cate and costly to manufacture. Pressure losses in these passages are 
high and must be known in order to establish pump and turbine design 
point data. Accurate prediction of these losses by analytical means is 
not possible except for flow passages with simple geometries. There- 
fore, it is advisable to conduct model tests with fluids at low pres- 
sure, preferably air, to determine pressure loss coefficients experi- 
mentally which can then be used to predict the pressure losses in the 
actual engine. 

The ARES (Advanced Rocket Engine Storable) propulsion unit is a 
compact dual -combust ion liquid rocket engine designed by the Aerojet- 
General Corporation, Sacramento, California. Professor M. H. Vavra of 
the Department of Aeronautics, Naval Postgraduate School, Monterey, 
California, has designed full-scale fluid dynamics models of this unit 
and has conducted preliminary tests of the oxidizer flow passages with 
air. The models have been made available to the writer by Aerojet- 
General and were used to carry out more extensive and more detailed 
testing with improved instrumentation. 
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The purpose of this thesis is threefold. First, the air test 
results are used to predict the pressure drops in the flow passages, 
flow distribution within the unit, and velocity distribution at the 
entrance to the combustion chamber for a liquid flow. Second, the 
effects of compressibility and Reynolds number on the test results are 
investigated, and the validity of the testing method is established. 
Third, changes in design are recommended based on the test results and 
flow visualizations. 

The author wishes to acknowledge the guidance provided by Profes- 
sor Vavra and extends sincere thanks to Mr. R. W. Savage for his assis- 
tance in performing experiments. Acknowledgment is due also to Mr. R. 
Beichel, Aerojet-General Corporation, for the loan of the models and 
his endeavors to apply model testing techniques to the design of ad- 
vanced rocket engines. 
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SECTION 2 



DESCRIPTION OF ARES UNIT AND TEST MODEL 

A three-dimensional schematic representation of the oxidizer 
passages of the ARES propulsion unit is shown in Fig. 1. The oxidizer 
is composed of a mixture of nitrogen tetroxide and hydrazine and it is 
used for regenerative cooling of the combustion chamber and the rocket 
nozzle before being mixed with the fuel. Oxidizer enters the outer 
annulus through a centrifugal pump. This annulus contains turning 
vanes, not shown in Fig. 1, to distribute the flow uniformly to the 
tubes that form the outer skin of the nozzle. Each tube is an inte- 
gral unit consisting of a down-flow and an up-flow passage..^ The en- 
trance and exit holes are in the same plane perpendicular to the nozzle 
axis; however, they are located at different radii in this plane. The 
tubes begin and end with an inner diameter of 0.3 inches. At the noz- 
zle throat they form an ellipse with major and minor axes of 0.31 in. 
and 0.19 in., respectively. This elliptical cross section gradually 
increases until at the bottom of the tubes the major and minor axes 
have values of 0.76 in. and 0.70 in. A scale model of one of these 
tubes is shown in Fig. 2. From the discharge of the tubes the oxidizer 
enters a second annular chamber concentric with that into which the 
pump discharges. In this annular passage the flow is turned to enter 
a series of 60 holes drilled radially in the inner shell of the second 
annulus. The oxidizer is then discharged into the primary combustion 
chamber through 180 injectors. A flow passage, not shown in Fig. 1, 
directs the combustion gases to a single-stage axial turbine which 
drives the pump previously mentioned. Details of this flow passage are 
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given by Vavra. Tests of two turbine designs have been performed, the 

2 3 

results of which are reported by Vavra and Messegee. 

The portion of the unit from the inlet of the outer annulus to the 
bottom of the rocket nozzle is called the oxidizer discharge. The por- 
tion including the remainder of the nozzle tubes to the injectors is 
called the oxidizer return. Secondary combustion takes place after the 
turbine discharge and ahead of the nozzle throat, 

A scale model of these flow passages has several modifications as 
far as shape and layout of the actual unit is concerned, but mostly re- 
tains flow passage dimensions. The model, shown assembled in Fig- 3, 
is constructed of plastic and aluminum and weighs 136 pounds. A sche- 
matic of the flow passages is shown in Fig. 4. The discharge and re- 
turn annuli are separated by 7 parts of varying flow areas and shapes. 
Flow enters the model through simulated pump vanes which impart a whirl 
component. The whirl is removed by turning vanes and the flow is dis- 
tributed to part 1, shown separately in Fig. 5. Part 1 has 32 equally 
spaced groups of 3 holes with a total flow area of 4.712 sq. in. Parts 
2 and 3 are shown respectively in Figs. 6 and 7. Each of these parts 
has 52 equally spaced single holes with a total flow area of 2.552 sq. 
in. Parts 4, 5, and 6 together with 2 orifice plates simulate the noz- 
zle tubes. An exploded view of this assembly is shown in Fig. 8. The 
first orifice plate, having holes of diameter d * 0.217 in. and dia- 
meter ratio d/D 3 ^ = 0.723 simulates the dpwtilfiosr ofufehfi noizlfe 
tubes. The second plate with d = 0.216 in. and |3 = 0.720 simulates 
the up- flow of the nozzle tubes. Part 7, s'hown in Fig. 9, has 32 
equally spaced triple holes with a total flow area of 5.574 sq. in. 

From part 7 the flow passes into the second annulus shown in Fig. 10, 
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where it is turned to enter the injectors. Each of the 60 injector 
tubes has 3 holes labeled A, B, C, respectively, from the outer radius 
inward. The tubes are labeled clockwise from 1 to 60 with the first 
tube just past the nine o'clock position when the injectors are viewed 
from the combustion chamber (see Fig. 36). 

The model can be divided into sub-assemblies for testing. Tests 
on the oxidizer discharge sub-assembly including the first annulus 
through part 5, excluding the second orifice plate, were made by Vavra. 
Preliminary tests of the return sub-assembly, part 6 through the in- 
jectors including the second orifice plate, were also made by Vavra. 
This paper presents more complete tests of the return sub-assembly, in 
addition to tests of the fully assembled unit. 
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SECTION 3 



TEST INSTALLATION AND INSTRUMENTATION 

Available as a source of air is a 200 H. P. centrifugal compres- 
sor which delivers about 4 pounds mass of air per second at a pressure 
ratio of about 2. Part of this flow passes through a cooler; the rest 
discharges through a butterfly valve to atmosphere. Figure 13 is a 
schematic of the flow measuring installation. The cooled air enters 
plenum 1 then passes through a 4 in. pipe discharging into plenum 2 
shown in Fig. 11. Flow measurement is accomplished by either a sharp 
edge orifice or a flow nozzle installed in the 4 in. pipe (Fig. 12). 

The design of the orifice and the flow nozzle is in accordance with VDI 
standards^, and the installation satisfies or exceeds the specified re- 
quirements of pipe lengths upstream and downstream of it. In addition 
a flow straightener is installed at the discharge of plenum 1. The 
orifice has a diameter d * 3.024 in. and a diameter ratio j?> = 0.7429. 
The nozzle has d m 3.022 in. and =■ 0.734. Initial tests conducted 
by Vavra used the sharp edge orifice to measure flow rate. Prior to all 
tests performed by the writer the flow nozzle was installed. 

Appendix A contains the derivation of the flow measuring equation 
for the nozzle. Similar information on the orifice is given by Kelly. ^ 
Plastic hoses connect pressure taps with the measuring instruments. 

The pressure difference P - P^ is measured by a Meriam Vernier Mano- 
meter filled with water. The gage pressure ahead of the nozzle is 
measured by a similar manometer filled with mercury. The measuring 
accuracy of these manometers is about + Q.Q02, :j.n.. The temperature t 
is measured with an accuracy of + 1°F. by an iron-cons tantan thermo- 
couple located 48.5 in. upstream of the nozzle. 
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Within plenum 2 there is arranged a bend and a 10 deg. conical dif- 
fuser with a series of screens at the discharge to minimize flow insta- 
bilities. A second thermocouple is situated within plenum 2. Plenum 
pressure is read on either a mercury filled manometer or a water filled 
manometer bank. Three such manometer banks were used. Bank 1 is a 
250 cm. 40-tube bank which can be read to an accuracy of + 0.1 cm. 

Bank 2 is a 175 cm. 60-tube bank with a reading accuracy of + 0.2 cm. 
Bank 3 is a 92 in. 20- tube board that can be read with an accuracy of 
+0.1 in . 

The model is attached to plenum 2 by means of a steel plate visi- 
ble in Fig. 11. Mounted on this plate is a wooden nozzle of circular 
cross section to produce uniform entrance conditions. Situated at the 
throat of the nozzle is a Kiel probe measuring total pressure P 

11 in 

For tests with the entire model assembly the conduit shown in Fig. 14 
is used to introduce air from opposite directions into the first annu- 
lus. Static pressure taps measuring P q are located 1.5 in. ahead of 

the model inleto In the same plane are located probes to measure P fc . 

o 

One such probe is pictured in Fig. 15 and shown schematically in 
Fig. 16. The design allows rotation of the probe into the flow, en- 
suring an accurate total pressure reading. 

Figure 4 shows the location of the measuring planes within the 
model. Pressure taps measure static or total pressure depending on 
their location. In each of the planes Z, B, C and E are arranged 
8 equally spaced static taps. Planes Y and G each have 8 equally 
spaced total taps. Planes A and F have 4 static and 4 total taps, 
spaced alternately. Plane H has 12 equally spaced static taps. Plane 
X contains no taps, but the pressure can be accurately predicted at 



21 



this location. Total and static pressures are measured at selected 
injector nozzles in plane I by the probes shown in Fig. 17. These 
probes, having a 15-deg. chamfer at the inlet, were designed to mea- 
sure total pressure of flows with yaw angles up to 20 deg. From 
Vavra's results it was suspected that the flow did not leave each in- 
jector axially. Therefore the new probe assembly was designed to per- 
mit quick adjustment to reach the condition depicted in Fig. 18. Total 
and static pressures at seven locations within the second annulus were 
obtained by inserting a removable probe which could be oriented into 
the direction of flow. All pressure taps were connected to the mano- 
meter bank by plastic tubing and measured against atmospheric or ple- 
num pressure depending on the test configuration. 

Testing of the return sub-assembly required a modification of the 
installation. An adapter fitted with screens, shown in Fig. 20, was 
mounted on plenum 2. Attached between the adapter and the sub- 
assembly was a ring of 52 inlet nozzles to introduce the flow uniformly 
to the orifice ring. Plane D, only used with this configuration, has 4 
pressure taps. Because flow is essentially stagnant within the adapter, 
it is assumed that =» . A view of this test installation, without 

D 'd 

the pressure taps connected, is shown in Fig. 21. 

The inlet pressure for both configurations was remarkably steady 
and, except where noted, the pressures within the model did not fluc- 
tuate more than + 0.5 cm. of water even at the highest flow rates. 
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SECTION 4 



THEORETICAL CONSIDERATIONS 

If a steady incompressible, adiabatic, viscous flow exists in a 
pipe, Bernouillis equation between two points on a streamline can be 
written as 

v 2 V 2 

( 144 )?! + e 2J-+ e z l ■ < 144 > P 2 + ? 2i“ + <e z 2 + \ j.j (1) 

where : 

2 

P * static pressure (lb/in ) 

£ ■ density (lb^/ft^) 

V =* velocity (ft/sec) 
g - 32.17 (ft/sec 2 ) 
z * height above reference (ft) 
h^ * head loss (lb/ft^) 

Since the first three terms on each side of equation (1) represent 

the total pressure, the head loss is equal to the total pressure loss 

AP in the pipe. 

t l-2 

Experience shows that this loss can be expressed by the function 

4P - f( p , u , V, X., d, e) (2) 

1-2 v 

2 

where: ® viscosity (lb:sec/ft ) 

Ji * characteristic length (ft) 
d ■ characteristic diameter (ft) 

e =» roughness parameter based on average variation of inside 
radius (ft) 

Buckingham's 7T theorem^ states that the number of independent 
dimensionless groups used to describe a phenomenon with n vari- 
ables is equal to n-r where r is the number of basic dimensions. For 
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the basic dimensions of force F, length L, and time T, a dimensional 
analysis of equation (2) gives 



AP. 



' 1-2 



f v 2 /g 




(3) 



Equation (3) has four independent parameters 



AP 



' 1-2 



e v 2 / g 



£ Vd 



£ 

d 




Euler number 



Reynolds number ■ Re 



relative roughness 
geometrical factor 



Usually equation (3) is expressed by 




where ^ is the friction loss coefficient that depends on the relative 
roughness and Re. In pipe flow with Re higher than about 10^ the 
coefficient remains essentially constant. In addition, the so- 
called momentum loss occurs if pipes have either sudden contractions 
or expansions. It is primarily due to eddying turbulence and has been 
found experimentally to satisfy the relation 



AP. 



“ 1-2 






(5) 



where \ is the momentum loss coefficient which depends primarily on 
flow geometry and remains practically independent of Re. 
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In a flow channel where frictional and momentum losses occur, the 



two loss coefficients may be combined to give 

- y $4 { <« 

where Y is a total loss coefficient which depends on Re. If Re is 
sufficiently high, the coefficient is a function of geometry 

alone. 

The analysis can be extended to systems of flow passages by con- 
sidering an equivalent pipe for which equation (6) is valid. If such 

z 

a system is modeled to scale, the factors Tj/* and — are contained in 
one loss coefficient ^ . For air as a testing fluid, care must be 
taken to maintain low Mach numbers M in the model as equation (6) holds 
for incompressible flows only. 

For air at low M and with the velocity at point 1 as a reference 



_ w 



e A i 



where: w = flow rate (lb /sec) 

m 



■ area (ft ) 



Equation (6) may then be written as 



AP. 



. 2 



= % ? _SL 

1-2 * J 1-2 A 2 

e A 1 



Using the equation of state for a perfect gas 



e - 



'1 (144) 



RT. 



(7) 



( 8 ) 



(9) 
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where R * 53.35 



ft-lb 

lb °R 
m 



T “total temperature ( R) 
1 

equation (8) becomes 



Ap. 



' 1-2 



-4 F 

2gA . L > 1-2 



RT\ 



w 



( 10 ) 



The so-called equivalent flow rate referred to total conditions at 
1 is defined by 




(U) 



Equation (10) then becomes 




1-2 




( 12 ) 



whore - f 1 . 2/2Ai 2 

* 

w will be defined as the equivalent flow rate referred to static 
conditions. Unless otherwise stated w refers to total conditions. 

For a liquid of specific weight ^ (lb/ft ) equation (6) becomes 

Li 



AP,. - % 
1-2 

Regrouping terms 






Ap. 



' 1-2 



1-2 



gA 2 ^L 



.2 

w 



Tl 



(13) 



(14) 
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Determination of k^ ^ by air tests of amodel makes it possible to 
predict pressure losses in a unit for particular fluids at a given flow 
rate* 



Prior to the application of this theory to the ARES model it was 
decided to construct an arbitrarily shaped flow passage which could be 
tested with air and water varying Re and M. This passage, which is 
shown in Fig. 22, has both friction and momentum losses. It was de- 
signed to operate in the same Re range for both fluids. Its inlet and 
discharge areas are equal to 4 sq. in. At the inlet is fitted a con- 
toured nozzle to minimize entrance losses into the first section of the 
pipe. Internally the model has a sudden expansion, a flow s traightener, 
and a sudden contraction. The minimum flow area of 2.23 sq. in. occurs 
in the flow s traightener . 

Air was obtained from the installation previously described. The 
model was attached to plenum 2, as shown in Fig. 23. Inlet total pres- 
sure is measured by the Kiel probe, and inlet temperature was taken as 

that measured in plenum 2. Discharge pressure P was measured by 

ais 

one of the probes shown in Fig. 17. Since conditions at the discharge 
were found to be nearly uniform, except in close proximity to the walls 
of the pipe, the average of five traverse points was taken as P 



t dis 



Static pressure P^^ was taken to be atmospheric. 



Values of k, , . , inlet to discharge pressure loss coefficient, 
m-dis 

were established for Re (based on d^) between 1.58 x 10"* and 3.67 x 

10 5 and are plotted in Fig. 24 with the curve giving the best fit 

through the data points. Variation of k^ n _^ g is from -0. 042 to 

+0.057 k . Higher values of Re could not be reached due to a lim- 
avg 

iting manometer height. The highest Mach number M, obtained at the 
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discharge as determined by the static to total pressure ratio, was 

O. 28. A slightly higher M was undoubtably present in the flow 
straightener . 

Figure 25 shows the test installation for the water tests of this 

model. At the 4 in. discharge of a fire hydrant were available about 

400 gal. l^O/min. at a pressure of 50 psia. After steady flow was 

reached, the discharge was deflected to fill a 50 gallon drum. By 

weighing the drum before and after each run and knowing the elapsed 

time, the flow rate was found. Inlet static pressure was measured 

in 

just after the entrance nozzle by a 50 psig. pressure gage.. 

5 

Values of k. were calculated for Re between 1.33 x 10 and 

m-dis 

4.88 x 10“* for 4 separate runs. The method of data reduction is given 
in Appendix B, and the results are shown in Fig. 24. Except for a few 
points there is excellent agreement between water data and the curve 
established by the air test. A maximum deviation of 10 per cent is 
noted while most points fall within 2 per cent of the curve obtained 
from the air tests. At both extremes of the Re spectrum the pressure 

P. tended to fluctuate, which accounts for some of the scatter in 
in 

these areas. What occurs at higher Re is hard to predict although it 
appears that k^ ^ is tending to approach some constant value. 
Because of limitations imposed by the method used to collect and weigh 
water, the higher Re range could not be tested. If a flow measuring 
orifice were installed upstream of the model inlet so that high flow 
rates could be measured accurately, data at high Re could be obtained. 

These tests prove that the use of equations (12) and (15) for the 
evaluation of the pressure loss coefficient for liquid flows from air 
tests is justified. 
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SECTION 5 



DESCRIPTION OF TESTS 

In addition to the two configurations of the ARES model and the 
flow channel mentioned previously, tests were made on a full scale 
model of a nozzle tube to determine a representative loss coefficient. 
Moreover, separate tests were made with a geometrically similar model 
of a tube section in the ARES model from plane W to plane X, simu- 
lating the nozzle tubes. These tests are described in the sequence 
they were made. Table I gives an index of the runs conducted. All raw 
data from these tests are on file with Professor Vavra at the Pro- 
pulsion Laboratory. 

5. 1 Oxidizer Return Sub-assembly 

In testing the oxidizer return sub-assembly, a total of six 
runs was carried out at Re between 0.868 x 10^ and 1.343 x 10^ referred 
to conditions just upstream of the flow measuring orifice. The maximum 
Re attainable was restricted by the height (250 cm) of manometer bank 1. 
As a result of Vavra 1 s preliminary tests it was determined that a new 
discharge total probe assembly should be designed. Runs Rl and R2 
were conducted without taking any measurements in discharge plane I. 

The new probe assembly shown in Figo 17 was built, calibrated, and 
used first during run R3 where alignment and ease of adjustment were 
investigated. More complete data in plane I were taken during runs 
R4, R5 and R6. At each selected peripheral location the probes were 
aligned individually with the injectors. Although the readings were 
steady and accurate, the procedure was time consuming. It was therefore 
decided to use data from every third injector tube only. 
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Pressure taps were made in the annulus by drilling a 1/16 in. 
hole partially tapped with a 10-32 thread. If the portable probe shown 
in Fig. 19 was not inserted, these holes were plugged by screws. 

Because of wear on the threads in the plastic, conditions in the annu- 
lus were determined only for runs R3 through R6. 

Pressures from plenum 2 downstream through plane H were measured 
on manometer bank 1 referenced to atmosphere. Pressure in plane I was 
measured on bank S^also referenced to atmosphere. Because of the 
large number of tubes and the slight pressure fluctuations noticed dur- 
ing testing, Polaroid pictures were taken of bank 1 for more accurate 
data recording. 

At locations 5, 6, 7 and 8 in plane H the pressures were observed 
to fluctuate considerably, especially at higher flow rates. The area 
of this unsteadiness is visible in Fig. 10. For runs R1 and R2 the 
fluctuation was dampened by clamping the pressure lines leading from 
these taps to the manometer bank. During runs R3 through R6 the clamps 
were removed and an average was taken between the maximum and minimum 
fluctuation. This method produced the same readings as with the clamps, 
but the fluctuations could now be observed. 

The general procedure for conducting each run was as follows: 

1. The micromanometers and thermocouples were zeroed and 
atmospheric pressure recorded. 

2. The desired pressure in plenum 2 was set and the flow was 
allowed to stabilize. 

3. Flow measuring data were recorded and a photograph of bank 1 
was taken. 

4. Plane I and annulus data were recorded. 

5. Part 3 was repeated. 



30 



Because of the interesting velocity distribution in the discharge 
plane and the large pressure fluctuations at locations in plane H, 
attempts were made to observe the flow in the annulus visually. Dif- 
ferent methods to introduce smoke had negative results. The first 
attempt used a 50 gallon drum within which a standard Navy MK-5 day 
signal flare was ignited. This particular flare produces about 15 
minutes of a dense, slightly toxic, white smoke. The drum was con- 
nected directly to plenum 2 and pressurized so that the smoke would go 
directly into the plenum. At the discharge of the model there was 
arranged a device to pipe the flow out of the building. It was found 
that too much smoke was produced, and the mixing in plenunl 2 was so com- 
plete that there was no streamline definition in the annulus. The 
flare produced enough heat to melt the rubber pressure seals in the 
drum, allowing toxic smoke to enter the room. For these reasons this 
attempt was abandoned in favor of introducing; smoke directly into the 
annulus. Cigars soaked in mineral oil and paraffin were ignited in a 
small container which was pressurizedo Lines from this container led 
to eight 1/16 in. diameter holes spaced around the entrance to the 
annulus. The smoke partially condensed in the lines and that part which 
entered the annulus diffused quite rapidly. Although streamlines were 
faintly discernible, they were of insufficient definition to observe the 
flow. 

This latter method left a large residue of tar and nicotine in the 
model. With the air passing through water was introduced into the 
annulus to clean out this residue. It became evident then that water 
would work as a visualization fluid. The set-up shown in Fig. 21 con- 
sists of a manifold connected to a plastic liquid soap bottle. Water 
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colored by food dye was forced into the annulus through the top four 
holes by squeezing the soap bottle. This method was very successful 
and produced a flow picture of photographic quality. Motion pictures 
of runs R3 through R6 were taken using a Kodak Cine-Special II, 16 mm 
camera. Both black and white and color film were used at 64 frames per 
second. Runs R3, R4 and R5 were photographed in black and white using 
TRI-EX, ASA 120 film and a red filter. KODACHROME II, ASA 25, was 
used on runs R5 and R6. 

5.2 Nozzle Tubes 

The nozzle tubes were tested in the configuration shown in 

Fig. 26. Water was supplied by a 1*0 in. hose from a faucet where the 

pressure was higher than 50 psig. A total of 22 data points during two 

4 4 

runs was taken between Re of 1.844 x 10 and 6.874 x 10 , referred to 
inlet diameter. 

Static pressure was measured at the tube inlet by a 50 psig pres- 
sure gage accurate to within + 0*1 psi. The flow discharged to atmo- 
sphere* For each data point the inlet pressure was set and the flow 
allowed to stabilize. Water was then collected in a five gallon bucket 
during a measured time. Data recorded were inlet static pressure, 
weight of water, and elapsed time. The maximum Re was obtained with 
the water supply fully open. 

5. 3 Geometrically Similar Model of Plane W to Plane X 

A single tube from plane W to plane X was scaled up by a fac- 
tor of 3.5 and tested with water by using the set-up shown in Figs. 27 
and 28. Water was supplied through a 1 in. hose to a 24 in. length of 
6 in. diameter pipe, at the bottom of which is installed a contoured 
nozzle to ensure smooth entry into the 1.05 in. diameter pipe section. 
For runs W1 and W2 , the orifices A and B had diameters of 0.761 and 
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0.757 in., respectively. Prior to run W3 new orifices, each having a 
diameter of 0.714 in., were installed. A constant water level in the 
reservoir was maintained at all times, and the flow rate was con- 
trolled by regulating the water supply and the discharge valve. The 

Reynolds number Re, referred to inlet diameter, could be changed from 
4 4 

0.415 x 10 to 4.078 x 10 during the tests. The runs were conducted 
in the same manner as those for the nozzle tubes. The pressures at 
the 6 locations indicated in Fig. 27 were read from a manometer board 
referenced to atmosphere. 

5.4 Discharge and Return Assembly 

Based on the information from these water tests two new ori- 
fice plates, with holes of 0.204 in. diameter, were installed in the 
ARES model. A total of five runs, R7 through Rll, was carried out at 
Re, referred to the 4 in. supply pipe, between 0.917 x 10"* and 1.32 x 
10^. Pressures from plenum 2 downstream through plane B were measured 
on bank 1 referenced to plenum 2. Bank 2, referenced to atmosphere, 
measured pressures from plane C through plane He Bank 3 was used again 
for plane I. The test procedures were similar to those used for the 
return sub-assembly. Two cameras were used to photograph manometer 
banks 1 and 2 simultaneously. No data was taken within the second 
annulus. Total pressure in plane Z was measured by the removable probe. 
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SECTION 6 



RESULTS AND DISCUSSION 

The oxidizer return sub-assembly was tested separately, mainly to 
investigate the velocity distribution in the discharge plane and the 
stability of the flow inside the return annulus. To a large extent 
these criteria determine the feasibility of the flow passage design. 
Assuming that fuel is pumped into the combustion chamber at plane I and 
distributed uniformly around the periphery, the importance of a uniform 
oxidizer distribution can be seen. The mixture of fuel and oxidizer 
will be rich or lean depending on whether the oxidizer flow is defi- 
cient or excessive. If the oxidizer discharge varies either with time 
or peripherally, and if it is not possible to meter the fuel accord- 
ingly, a non-uniform combustion will result. This condition is unde- 
sirable for two reasons. First, hot spots will result in the combus- 
tion chamber inducing thermal stresses in the combustion liner. Second, 
and more important, is the effect on turbine inlet conditions. The 
combustion chamber length is short, and non-uniform temperatures would 
have insufficient time to adjust prior to reaching the turbine inlet. 

A turbine blade then would pass through a non-uniform temperature field. 
An off-design condition would result and could have a detrimental effect 
on turbine performance. A major portion of the regenerative cooling in 
the ARES system occurs within the return annulus. The convective heat 
transfer in this annulus is totally dependent on the oxidizer flow 
pattern. 

Since the return sub-assembly differs from the full model only 
upstream of plane E, it seems reasonable to assume that flow downstream 
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of this point will represent the conditions that will exist in the 

fully assembled model. Actually it was found that inlet conditions at 

plane E were essentially the same in both test configurations. 

Data reduction for runs Rl through R6 is contained in Appendix C, 

* 

part 2. The equivalent flow rate w^ referred to conditions at the 

flow measuring nozzle, and the values of Re just upstream of this point 

are listed in Table II. Appendix A gives the derivation of the flow 

measuring equations. Actual flow rate w varies from 0.283 lb /sec in 

m 

run R3 to 0.449 lb /sec in run R6. Also listed in Table II are average 
m 

total pressures and equivalent flow rates at intermediate planes D, F 
and G found according to equation (C2) . Losses between planes will be 
referred to one of these average inlet conditions. 

It is important that the flow is distributed evenly around the 
periphery as it enters the model. In this case the procedure of aver- 
aging conditions at each plane can be justified. From water tests con- 
4 

ducted by Vavra , on a geometrically similar model of a tube from plane 
D to plane X, an orifice discharge coefficient K was obtained. Vavra 1 s 
model scaled the orifice size by a factor of 3.5 from 0.216 in. to 
0.757 in. This discharge coefficient, analagous to oC f r of equation 
(A3) , takes into account some of the difference between measured and 
actual flow rate. Figure 29 shows K plotted as a function of Re 
referred to tube inlet conditions. Values of K in Fig. 29 appear some- 

g 

what higher than those listed by ASME ..for a sharp edge orifice with 
^ = 0.7209. This is quite reasonable when the different inlet con- 

ditions are considered. The ASME data are dependent on a minimum 
length of straight pipe upstream of the orifice. For Vavra 1 s tests an 
inlet nozzle was present only ahead of the orifice. In later tests of 
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the same orifice with some length of pipe upstream, values of K were 

obtained which agreed more closely to those listed by the ASME. 

The flow at plane D was found to be stagnated, and the pressure 

P = P was uniform. At each measuring point in plane E an equiva- 
to * 

lent flow rate w^ was found according to equations (C5) and (C6) . 
**■ 

Assuming that w^ existed through each of the 52 orifices the peri- 
i 

pheral flow distribution can be expressed by a factor F^, where 



(52) w 



F. = 

1 



D. 

l 



w^ 



(15) 



The total equivalent flow rate w^ is given in Table II. Tables III 
and IV list values of F^ in addition to intermediate quantities used 
in the calculation. Figure31 is a plot of F^ as a function of peri- 
pheral location for Rl through R6. Looking downstream the angle 
is measured counterclockwise with 0=0 deg. corresponding to the 
three o'clock position. In each instance the flow distribution is 
almost identical. The maximum deviation of this average F^ from 
F^ =* 1.0 occurs in run R3 where the lowest Re was encountered. It was 
observed that most of the questionable data obtained during these tests 
occurred at the lower Re where the measuring accuracy is lower. Gen- 
erally each curve has two peaks and two valleys indicating a deficiency 
of flow at 0 = 110 deg. and 300 deg 0j and an excess at 0=0 deg. and 
225 deg. This is a relatively symmetric condition with the maxima 
appearing about 180 deg. apart. Because of this somewhat regular 
occurrence the values of F^ at © = 90 deg. for runs R3 through R6 are 
questionable. Not be be overlooked, however, is the fact that the 
maximum deviation from F^ ■ 1.0 is only 4 per cent. This fact tends to 
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establish the accuracy of Fig. 29 and to substantiate that flow enters 
the oxidizer return sub-assembly uniformly. 

Figure 30, also plotted from Vavra's water tests, shows ^ , 

U"A 

the loss in total pressure from plane D to plane X divided by the static 
pressure change across the orifice, as a function of Re. The total pres- 
sure P at plane X is found according to equation (C8) . With this 



pressure and the inlet conditions at plane D, the coefficient ^ was 
calculated for each peripheral location. An average value for 

each run was found which is shown in Tables III and IV. It should be 
noted that at any peripheral location the value of k^ differs from 
k^ X by not more than five per cent. 

If, as was hypothesized, the losses in this section are indepen- 
dent of Re, then k^ should not vary from point to point. This may be 

shown by referring to equation (12) and recognizing that both AP 

* 2 _ 
and Wp depend essentially on The distribution of 

follows the same pattern as that established for F^. This is shown in 

* 



Fig. 31 for Rl and R2 and holds for the other runs also. Since w is 

l 

proportional to Re in accordance with equation (C6) , the loss coeffi- 
cient is also proportional to Re. 

Experimentally obtained values of the friction factor proportional 

to k as a function of Re and relative roughness are shown in Fig. 32 

for two different conditions of flow in pipes. These curves were taken 
o' 

from Schlidhting- For the artificially roughened pipe of Fig. 32(a) 
the friction factor increases with Re to approach a constant value, 
whereas for the commercial pipe of Fig. 32(b) the factor decreases with 
Re toward a constant value. In the former case the pipes tested were 
lined with protuberances which penetrated completely through the 
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laminar sub-layer of the boundary layer, A large part of the friction 
is then due to a form drag which can be thought of as a momentum loss. 
In the second case skin-fric tion alone accounts for the friction fac- 
tor. As relative roughness is increased in both cases the friction 
factor reaches a constant value at a lower Re. Realizing that in the 
present model test each tube from plane D to plane X is machined to at 
least the equivalent of commercial roughness and that k^ ^ increases 
slightly with Re, it can be concluded that the momentum loss created by 
the orifice predominates. Tables III and IV show values of Re corres- 
ponding to the general transition range of Fig. 32. The loss coeffi- 
cient k^ would become constant if either more momentum loss were in- 
troduced or if Re were increased. 

Incompressible flow is a basic assumption of the theory proposed 
for this analysis. It seems reasonable to assume that the highest Mach 
number M in the model exists at the minimum flow area which occurs at 
the vena-contracta downstream of the orifice plate. At plane E an 



through (C13). Also calculated were average values of P £ , P and M^. 

E 

Table V lists these values for runs Rl through R6. Applying the values 



flow area varying from 2.081 sq. in. to 2.161 sq. in. The next 
largest area in the model is the 3,138 sq. in. area at plane I, where 



the highest M in the models occurs at plane E, reaching a maximum value 
of 0.34 in run R6. Hence it is permissible to assume that the flow is 
essentially incompressible. 



average area restriction factor 




found using equations (C9) 



of ^ to a flow area at plane E of 3.676 sq. in. gives an effective 
E 



an average M was calculated and is listed in Table V. For all runs 
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Table V also lists values of average loss coefficients from planes 



D to F) F to G, G to I and X to F. The latter loss is given by 



X-F 




(16) 



From Fig. 30 are obtained values of ^ which can be used to pre- 
dict , and hence whereas the values are calculated from 

X 

test data. The magnitude of k is small because of the small dis- 

X* F 

tance between the two planes. It can be said, therefore, that the 
method used to establish Fig. 30 is a valid one. 

The peripheral distribution of is shown in Table VI. 

Except for a few isolated values, is contained within 2 per cent 

of Referring to the previous discussion concerning the rela- 

tionship between equivalent flow rate and loss coefficient, it can be 
assumed that the flow is distributed uniformly at plane F. 

The situation is different in plane G. Figure 33 shows the per- 
ipheral distribution of k /k . Since the pressure tap at 0 3 0 deg. 

F“G F“G 

was damaged, a curve was faired in to give an estimated value at this 
location. Except for runs R3, at low Re, a definite pattern can be 
observed. A minimum loss coefficient occurs at 0 s 180 deg. This in- 
dicates an irregularity of the flow at this peripheral location. This 
situation can be explained from the geometry of the flow passages. 

Figure 34 shows the intersection of the two cylinders with the posi- 
tions of Q in plane G indicated. At © - 90 deg. and 270 deg. the 
flow has an easier path to follow, therefore one would expect a higher 
flow rate at these points and hence a higher value of loss coefficient. 
However, at & = 0 deg. the flow enters the return annulus on top of 
the holes in plane H. It will be shown than the primary injectors at 
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this location pass a higher percentage of the discharge flow. It is 
entirely reasonable then that the flow is deficient as it enters the 
annulus from Q = 135 deg* to Q = 225 deg. 

Table VII lists values of loss coefficient k for each location 

G-H 

in plane H and k for runs Rl through R6. Also shown in Table VII 
G”H 

are values of k^ ^ for runs R3 through R6, found by 
^H-I = ^G-I " ^G-H 

Shown in Fig. 35 are values of average loss coefficients between 
different stations of the model as a function of Re. The increase with 
Re of k^ ^ has already been explained, however, the maximum and minimum 
values are only 92 and 104 per cent of the average value, respectively. 
The highest value which k^ ^ approaches at increased Re is probably 
around 0.125. If loss coefficients in the low Reynolds number range 
are not .considered, v th6 average losses from plane X throughout the re- 
mainder of the model are remarkably constant. 

Velocity distribution in the discharge plane is found according 
to equations (C15) and (C16). Figure 36 shows a drawing of plane I 
from the discharge side and the relative position of the pressure taps 
in plane H. Table VIII lists values of V^/V for each of the surveyed 
locations for runs R3 through R6. An average of positions A, B and C 
was calculated at each location 1 through 57. This average velocity 
distribution is plotted in Figs. 37 through 40. The velocity distribu- 
tion and hence the flow distribution is essentially independent of Re. 
The excessive flow from locations 54 through 10 was anticipated from 
the conditions existing at 0=0 deg. in plane G* The erratic con- 
dition from locations 24 through 36 will be explained later. It is not 
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understood why the flow is deficient around location 48. Data at the 
seven locations within the annulus obtained from the removable probe 
were entirely unreliable. The model had been dropped causing damage to 
the plastic threads and making it virtually impossible to prevent 
leakage. 

The flow within the return annulus can best be described by re- 
ferring to the results of the flow visualization efforts previously de- 
scribed. It has been shown that the flow in plane G is symmetrical 
about the 0 = 0 deg. and 0 = 180 deg. axis. Therefore, the annulus 
is naturally split into a top and bottom half in each of which the same 
flow pattern should exist. Water introduced in the lower half had a 
tendency to pool at the lowest point because of gravity effects, vhereas 
all the water introduced in the top half was carried out of the model. 
For these reasons and because of problems associated with photographing 
the flow conditions, it was decided to investigate only the upper half 
of the annulus. Figure 41 is a drawing of that half annulus, developed 

into a plane, showing the locations where water was introduced. In 

4 

addition, data obtained by Vavra are shown that give flow directions 
at locations within the annulus. These flow directions are represented 
by arrows in Fig. 41. Figure 42(a) through (h) and Fig. 43(a) through 
(h) are photographs taken of runs R4 and R5, respectively. The upper 
half of the annulus is shown as the upper 3/4 of each picture. These 
pictures were made from single frames of 16mm movie film taken at 64 
frames per second. Each sequence spans one second of operation with a 
time between each pattern of 0.125 sec. Figure 42 is typical of the 
flow pattern at Re 53 1.028 x 10^ and lower whereas Fig. 43 is repre- 
sentative of the conditions at Re = 1.177 x 10^ and higher. Several 
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similarities exist at all Re. Most apparent is the extreme instability 
in the region of plane H where the flows from the top and bottom halves 
of the annulus mix. This region is exactly that in which the pressure 
taps H5, H6, H7 and H8 are located. 

It should be noted that this is the only region in the model where 
the pressures were unsteady. Several vortices are present in this gen- 
eral area. The pattern is apparently being triggered by a huge separa- 
tion from the turning vane near H5. Although not entirely obvious from 
the sequence of pictures, it was observed that these vortices went 
through a periodic variation in both intensity and location. At best, 

it can be said that the flow in this area is thoroughly confused and 

that this situation is most likely the cause of the cyclic velocity dis- 
tribution in plane I which is evident from Figs. 37 through 40. It has 
been noted also that a deficiency of flow exists at the far right of 
the annulus. This portion of the annulus is supplied by the 9 = 180 
deg. location of plane G where the deficiency of flow has already been 
observed. The most striking difference between Fig. 42 and Fig. 43 is 

the appearance in the latter of a huge vortex in the upper right hand 

portion of the annulus. It is believed that this vortex is caused by a 
flow separation from the turning vane immediately upstream of it. A 
100 ft. roll of color film taken of runs R5 and R6 was left intact, and 
is on file with Professor Vavra. 

From the above discussion and the flow visualization studies, it 
is possible to make a qualitative design evaluation of this portion of 
the ARES assembly. The test data at the higher Re will best simulate 
actual operation, although this distinction appears necessary only when 
considering conditions within the annulus. Heat transfer will most 
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certainly be hampered where vortices exist, since the same flow parti- 
cles will tend to circulate in one spot instead of passing through the 
model. Hot spots could result in these regions. The deficiency of 
flow in the region farthest from plane H would have the same effect on 
heat transfer, only to a lesser extent. Elimination of these flow in- 
stabilities, in addition to improving heat transfer characteristics, 
would serve to reduce the losses from plane G to H and to distribute 
the discharge velocity more uniformly. One method of accomplishing 
this would be to redesign the turning vanes, making their shape conform 
more to airfoil cross sections cambered in the direction in which the 
flow is to be turned. For stress considerations, the contact surface 
area with the outer annulus would have to be maintained. Another 
scheme would be to replace the relatively large number of turning vanes 
with a few long flow dividers. It is felt that a large portion of the 
instability occurs as a result of the mixing of flows from each half of 
the annulus, and for this reason the arranging of a flow splitter is 
recommended. 

The loss coefficients calculated during the tests of the return 
sub-assembly were duplicated when the full model was tested and will be 
discussed later in this report. Before proceeding with the tests of 
the entire unit, it was decided to determine how closely the pressure 
drop in the model, from plane W to plane X, simulated that in the 
actual nozzle tubes. Results of tests on the nozzle tube are listed in 
Table IX. The method of data reduction is shown in Appendix C, part 3. 
Loss coefficient ^ ^ is plotted as a function of Re in Fig. 44. 

Fluctuation of the pressure gage was as much as + 0.5 psig. At Tofr 
values of Re, where pressures were of the order of 3.5 psig., this 
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fluctuation could cause the computed to vary between 11.1 and 8.5. 

Therefore, the dashed portion of the curve of Fig. 44 is questionable. 
The decrease of ^ with Re is quite evident, and by applying the 
rationale previously described, it seems that skin friction accounts 
for the major portion of the loss. Because of the scale chosen for ^ 
in Fig. 44 its variation appears magnified. Actually ^ only deviates 
4* 6 per cent from an average value of 10.3. The range of Re encountered 
falls well within the transition range for pipes as shown in Fig. 32. 

Data from tests on the geometrically similar section of plane W to 
plane X were reduced by the method outlined in Appendix C, part 4. 

These tests were conducted similarly to those run by Vavra for which the 
validity has been established. Table X lists values of ^ ^ ^ ^ ^ 

and the orifice discharge coefficients, in addition to other quantities 
of interest. Points 1 and 6 refer to planes W and X respectively. 

Figure 45 shows f and ? 1-6 £or runs W 1 and W 2 , plotted as a func- 
tion of Re. The loss coefficient ^ 1 6 by almost a factor of 

2 from ^ of the nozzle tubes. This condition is unsatisfactory 
because the largest loss in the entire assembly occurs in these tubes 
and the water model with orifice plates of d 3 0.757 in. and 0.761 in. 
gave half of this loss only. It is also interesting to compare the 
values of K (0.757) with those obtained by Vavra. The values differed 
on an average by almost 6 to 7 per cent. The difference in the dis- 
charge coefficient is due to the influence of the inlet conditions. 

Table XI also lists the results of run W3 with orifice diameters 
of 0.714 in., whereas Fig. 46 shows ^15 anc * i" 1 6 as ^ unctions of 
Re. It appears as if this orifice diameter more nearly simulates con- 
ditions in the nozzle tube. From the test of the fully assembled model, 
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equations (C32) to (C36) . Figure 47 shows these values together with 
the curve obtained from the water tests of the nozzle tubes. The dif- 
ference in inlet conditions to plane W when the section is installed in 



Apparently though the orifice size of 0.714 in. does the best job of 
simulation. 

Tests of the complete assembly of Fig. 4 were conducted mainly to 
establish values of loss coefficients without making pressure or veloc- 
ity surveys. The methods of data reduction are similar to those used 
for the return sub-assembly and are contained in Appendix C, part 5. 
Table XII gives the data for the flow rate determined, together with 
average total pressures and equivalent flow rates in planes Z, A, F 
and G. Data at plane Y indicated a lower total pressure than that 
occurring in plane A. This impossible situation must be due to inac- 
curate measurement in plane Y. These data were therefore disregarded. 
Average loss coefficients from planes Z to A, A to F, F to G, G to H, 

H to I and W to X are listed in Table XII for runs R7 through Rll. 

The range of Re encountered in the 4 inch supply pipe, and in the 
model at plane E, was the same as for the return sub-assembly alone ^as 
was the value of discharge Mach number M^o Arranging the discharge 
sub-assembly ahead of the return sub-assembly does not disturb the 
peripheral velocity distribution that was observed in the return sub- 
assembly alone. Table VI lists values of k * 7k „ which differ 

A - F A - F 

from a value of 1.0 by a maximum of 0.005. The distribution of loss 
coefficients from plane F to plane G is shown in Fig. 48. It is almost 
identical to that from runs Rl through R6. Velocity distribution in 



the model could account for the differences between 
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plane I was also found to agree precisely with that determined pre- 
viously. From the foregoing it is reasonable to conclude that the dis- 
cussion of results for the return sub-assembly is applicable to the 
entire model. 

Figure 49 shows the average loss coefficients as a function of 
Re. The scatter of ^ occurs because the total conditions at plane 
A were measured at only one location, hence a true average could not be 
obtained. However, this condition is not serious because this loss is 
only a small percentage of the total loss. The total loss coefficient 
j. is obtained by adding the coefficients between intermediate planes. 
If only the conditions in planes Z and I are considered, a value ^ * 
is obtained which differs from k^ As Re is increased k^ ^./k^ j* 

goes from 0.974 to 0.948. Evidently there exists a slight effect of 
compressibility. The difference in the two loss coefficients may be 
shown mathematically. Assume, for example, that the loss coefficient 
between planes Z and F is to be determined. The overall coefficient 
may be expressed as 



Z-F 



_* 2 - 

p. 



w, 



( p t - p t ) +(p t - p t ) ( 

Z A A F ) 



(18) 



Using the relation 



— * — — * — 

W z p t = W A p t 
L \ A C A 



(19) 



equation (17) becomes 



^-F 3 ^-A + P_ k A-F 



( 20 ) 
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Since total pressure must decrease from plane to plane, the pressure 



ratio P /P must be larger than unity. Hence, equation (20) shows 
_ z f A 

that will always be larger than The total loss coefficient 

from plane Z to plane I can be expressed by 



— > 



Vi " Va + k A-F + k F-G + k G-H + Vi 



( 21 ) 



where 



j 



A-F 




(etc.) 



Equation (21) was used to compute k^ for runs R6 through Rll. The 
values were found to be identical to those obtained from 




( 22 ) 



Consider the section from plane G to plane H only. Figures 35 

and 49 show that the loss coefficient as a function of Re, if based 

only on inlet conditions at plane G, is the same for both the return 

sub-assembly and the fully assembled model. If the method of equation 

(21) were applied, it can be seen that k 1 would increase as more 

G-H 

sections were added to the model upstream, because for the same inlet 

* 

conditions at plane G, namely, P and w , the pressure ratio P /P 

C G G C Z fc G 

would increase. This method actually refers all loss coefficients 

back to the initial inlet conditions and assumes incompressible flow 
from, say, plane Z to plane H. Using k to represent the loss, 
effectively cancels any errors introduced by the incompressible assump- 
tion for flows upstream. Only in that portioncof the model from 
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G to H is incompressible flow assumed. For this reason, it was decided 
to use the above-mentioned simple method for computing total loss 
coefficients. As an extension of this reasoning, the writer believes 
that more accurate estimates of loss coefficient could be obtained by 
adding more intermediate measuring planes. It is questionable whether 
this complication is justified since the maximum difference between 

y 

and is only 5.2 per cent. 

From Fig. 49 the coefficient ^ appears to approach a value of 

about 0.62 of which 66 per cent is from plane W to plane X. As an 

example it will be assumed that the flow rate of the oxidizer w is 

100 lb /sec and that its specific weight is 75 lb/cu. ft. The 
m L 

total pressure drop Ap in the oxidizer passage from the pump dis- 

Z-I 

charge to the entrance into the primary combustion chamber is then from 
equation (14) 




AP 



Z-I 



where: 



• Vi 

Vi " 



W 2 (144) 
g V L 



(lb/in 2 ) 



0.62 in 
32.17 ft/sec 2 



or 



(23) 



AP - °-62 (100) 2 (144) , 370 lb/ln 2 

Z-I (32.17) (75) 



(24) 
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SECTION 7 



CONCLUSIONS 

The average loss coefficient from a position just downstream of 
the oxidizer pump to the discharge of the primary injectors is 0.62. 

This experimentally determined coefficient should remain nearly con- 
stant near the specified operating range since the anticipated Reynolds 
number is sufficiently high. The model section from plane W to plane X 
with 2 orifice plates, having holes of d * 0.204 in., represents the 
actual nozzle tubes. The pressure loss coefficient found in the air 
test of this model section agreed closely with the coefficient deter- 
mined from the water test of the tubes. The peripheral distribution 
of flow within the model is very nearly uniform until it enters the 
annulus. Here the flow becomes deficient in a 45 deg. sector centered 
on 6 ■ 180 deg. in plane G. This causes a deficiency of flow in that 
portion of the return annulus farthest from plane I. 

Because of compressibility effects a difference exists between 
loss coefficients determined from inlet and discharge conditions only, 
and those obtained by a summation of the intermediate loss coefficients 
between adjacent planes. The flow channel used to check the basic 
theory by water and air tests was neither designed nor instrumented with 
this difference in mind. It would be of interest to construct a more 
complicated system of flow passages for this purpose and to investigate 
the effect of compressibility more fully. It is believed by the writer 
that more intermediate planes where exact conditions can be established 
will lead to more accurate estimations of loss coefficient. 

Due to the fact that the major portion of the pressure drop is 
caused by momentum losses, it would be inconsistent to apply standard 
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theories that relate heat transfer to pressure drop with the objective 
of evaluating the heat transfer coefficients at different locations in 
the flow passages. However, the vortices within the return annulus, as 
evidenced by the photographs of the existing flow pattern, will cer- 
tainly be detrimental to heat transfer in this region. These vortices 
could be largely eliminated by changing the design of the turning vanes 
and installing a flow divider where mixing now occurs between flows 
from the top and bottonrof the annulus. 

Model tests of this type provide a natural bridge between the 
completion of a design and the building of a prototype. It is estima- 
ted that a cost factor of 20 to 1 exists between the actual hardware of 
the ARES unit and its corresponding flow model. Prediction of the 
pressure loss coefficients of complicated passages analytically is, at 
present, an impossible task. Except for the slight effects of Reynolds 
number, the experimental method of predicting pressure loss coefficients 
presented in this paper is both valid and accurate. 
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TABLE I 





INDEX OF RUNS 


RUN NUMBER 


TEST DESCRIPTION 


RM1 


Air test of flow device 


RMWl, RMW2, 
RMW3, RMW4 


Water test of flow device 


Rl, R2, R3, 
R4, R5, R6 


Air test of oxidizer return sub-assembly 


Wl, W2 


Water test of geometrically similar model 
plane W to plane X^orifice plates of 
d ® 0.761 in. and d * 0.757 in. 


W3 


Same as Wl and W2 with orifice plates both 
of d = 0.714 in. 


WV1, WV2, WV3 


Water tests of geometrically similar model 
plane W to plane X 


R7, R8, R9 , RIO, 


Rll Air tests of full ARES assembly 
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TABLE II 

OXIDIZER RETURN SUB-ASSEMBLY, INLET ORIFICE 
DATA AND EQUIVALENT FLOW RATES 



RUN 

AVERAGE OF ALL 
DATA POINTS 


Rl 


R2 


R3 


R4 


R5 


R6 


<1Wi " 2) 


14.768 


14.768 


14.778 


14.778 


14.792 


14.734 


P u (lb/ in 2 ) 


17.758 


17.550 


16.094 


16.746 


17.437 


18.147 


A h Q (in H 2 0) 


1.998 


1.887 


0.971 


1.380 


1.776 


2.221 


T u (°R) 


539.0 


539.7 


540.5 


541.0 


552.0 


551.5 


T P len <° R > 


534.7 


535.3 


536.0 


542.0 


550.5 


549.5 


Re(1 ° = 4.07" 


1.307 


1.262 


0.868 


1.028 


1.177 


1.343 


★ 2 
W u (in ) 


.71785 


.70176 


.52655 


.61562 


.68344 


.74901 


P D (lb/in 2 ) 


17.708 


17.490 


16.058 


16.699 


17.367 


18.061 


P t (lb/in 2 ) 

F 


16.510 


16.385 


15.546 


15.926 


16.323 


16.704 


P (lb/in 2 ) 

G 


16.227 


16.122 


15.427 


15.738 


16.073 


16.383 


★ 2 
w D (in ) 


.71732 


.70137 


.52558 


.61799 


.68534 


.75132 


— * 2 
w p (in ) 


.76914 


.74870 


.54291 


. 64797 


.72918 


.81237 


— ★ 2 
» G (1- ) 


.78254 


. 76088 


.54708 


.65574 


.74051 


.82829 
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TABLE III 

FLOW DISTRIBUTION THROUGH ORIFICE PLATE, 
LOSS COEFFICIENT PLANE D TO PLANE X 

RUNS Rl, R2, R3 


Location 


1 


2 


3 


4 


5 


6 


7 


8 


Rl w * = .71732, K - .79, P " P. - .50 A P. = 1161.4 

u t D D-E 


p -p,, 

D E 

(Cm HjO) 


163.10 


159.40 


155.05 


151.80 


159.50 


163.10 


159.10 


158.45 


Re (I0' 4 ) 
d = 0.3 


3.503 


3.468 


3.426 


3.391 


3.468 


3.503 


3.466 


3.458 


* 

V (in 2 ) 


.01417 


.01403 


.01386 


.01371 


.01403 


.01417 


. 01402 


.01399 


F i 


1.027 


1.016 


1.005 


0.994 


1.005 


1.027 


1.016 


1.014 


^D-X ^ in ' ^ 


.1273 


.1244 


.1210 

avg. 


.1185 

.1239 


.1245 


.1273 


.1242 


.1237 


R2 w D * * .70137, K = .79, P t - P D - .50 A P = 1152.13 

X 


p -p„ 

D E 

(Cm H 2 0) 


151.25 


147.80 


143.85 


141.70 


147.65 


151.60 


147.60 


147.33 


Re(IO" 4 ) 
d = 0.3 


3.363 


3.328 


3.286 


3.264 


3.326 


3.366 


3.326 


3.324 


* 

\ (in 2 ) 


.01377 


.01363 


.01346 


.01337 


.01362 


.01378 


.01362 


.01361 


F i 


1.021 


1.010 


0.998 


0.991 


1.009 


1.022 


1.009 


1.009 


k D-X (in;4) 


.1250 


.1221 


.1188 

avg. 


.1171 

.1217 


.1220 


.1252 


.1219 


.1217 


R3 w * = .52558, K - .795, P„ - P_ - .505 A P. _ = 1095.29 

u t^ u D“E 


p -p 
D E 

(Cm HgO) 


69.70 


65.50 


68.95 


64.75 


66.85 


71.45 


67.55 


67.35 


Re (I0' 4 ) 
d =0.3 


2.249 


2.187 


2.238 


2.173 


2.207 


2.276 


2.218 


2.214 


w * 2 

W D. (in Z ) 


.01004 


.00976 


.00999 


.00970 


.00985 


.01016 


.00990 


.00988 


F i 


0.994 


0.965 


0.989 


0.960 


0.975 


1.006 


0.979 


0.978 


k D-X ,i "' 4) 


.1128 


.1060 


.1116 
av S- 


.1048 

.1096 


.1082 


.1156 


.1093 


.1090 
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TABLE IV 

FLOW DISTRIBUTION THROUGH ORIFICE PLATE, 
LOSS COEFFICIENT PLANE D TO PLANE X 

RUNS R4 , R5 , R6 


Location 


1 


2 


3 


4 


5 


6 


7 


8 


R4 w * = .61799, K = • 79 , P » P n ■ .508 A p » 1122.90 

D u u-fc 


P -P 
D E 

(Cm H 2 O) 


104.70 


99.50 


102.40 


96.60 


101.60 


106.70 


101.50 


100.50 


Re(l0' 4 ) 
d=0 3 


2.752 


2.687 


2.724 


2.650 


2.711 


2.780 


2.711 


2.687 


* 0 

D. (m ) 


.01188 


.01160 


.01176 


.01144 


.01170 


.01200 


.01170 


.01160 


F. 

1 


1.000 


0.976 


0.989 


0.962 


0.986 


1.008 


0.985 


0.981 


W"' 4 ’ 


.1185 


.1127 


.1160 

avg. 


.1094 

.1151 


.1150 


.1208 


.1149 


.1138 


R5 w * = .68534, K = .79, P. = P_ - . 50 A P^ _ = 1152.88 

D t D D-E 

A. 


p -p„ 

D E 

(Cm Hp) 


140.95 


135.05 


137.25 


132.05 


137.95 


142.95 


136.45 


135.95 


Re (I0' 4 ) 
d =0.3 


3.200 


3.139 


3.158 


3.106 


3.168 


3.220 


3.146 


3.146 


V (in 2 ) 


.01339 


.01313 


.0321 


.01299 


.01325 


.01347 


.01316 


.01316 


F. 

1 


1.016 


0.996 


1.003 


0.986 


1.005 


1.022 


0.999 


0.999 


k D-X (in ' 4) 


. 1228 


.1177 


.1196 


.1151 


.1202 


.1246 


.1189 


.1185 


R6 w * = .75132, K = . 79 , P = P - .50 A P. = 1181.06 

U t D D-E 

A 


VrS H 2 0) 


182.90 


176.05 


176.90 


171.35 


180.90 


185.60 


177.90 


177.05 


Re (I0‘ 4 ) 
d=0.3 


3.681 


3.618 


3.626 


3.574 


3.661 


3. 706 


3.606 


3.628 


* 

“d, (In 2 ) 


.01479 


.01454 


.01457 


.01436 


.01471 


.01489 


.01461 


.01458 


F i 


1.024 


1.006 


1.008 


0.994 


1.019 


1.030 


1.011 


1.009 




.1275 


.1228 


.1233 

av R- 


.1195 

.1245 


.1261 


.1294 


.1240 


.1234 
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TABLE V 

AREA RESTRICTION PLANE E, MACH # PLANES E, F & I, 

LOSS COEFFICIENTS PLANES D TO F, X TO F, F TO G & G TO I 


RUN 


Rl 


R2 


R3 


R4 


R5 


R6 


p t 

E (cm H 2 0) 


1169.33 


1159.50 


1096.98 


1125.75 


1155.63 


1184.63 


P E (cm H 2 0) 


1086. A3 


1082.89 


1061.75 


1072.87 


1084.22 


1091.77 


*E 


.32 


.31 


.21 


.26 


.30 


.34 


\ 

t d 


.979 


.981 


.991 


.986 


.982 


.977 


— * 

W E 

ST 


.81367 


.78919 


.55660 


. 67181 


.76516 


.86410 


<£ E 


.567 


.571 


.588 


.582 


.574 


.570 


P F (cm H 2 0) 


1148.37 


1140.39 


1087.82 


111.81 


1137.18 


1161.11 


P t 

F (cm H 2 0) 


1161.24 


1152.46 


1093.46 


1120.18 


1148. 11 


1174.90 


*F 


.12 


.12 


.08 


.10 


.12 


.13 


Vp (in' 4 ) 


.1309 


.1285 


.1166 


.1212 


.1280 


.1331 


Vf <in ’ 4) 


.0070 


.0068 


.0015 


.0061 


.0083 


.0086 


Vg (in " 4) 


.0293 


.0288 


.0261 


.0284 


.0290 


.0295 


P I (cm H 2 0) 






408.34 


407.81 


407. 71 


405.50 


I (cm H 2 0) 






420.92 


426.43 


432.78 


437.57 


“l 






.21 


.25 


.29 


. 33 


^G-I (in" A ) 






.0492 


.0501 


.0504 


.0517 
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TABLE m 



PERIPHERAL DISTRIBUTION OF LOSS COEFFICIENT 
PLANES A TO F AND PLANES D TO F 



K d -f 

K 0 - f 


RUN 


Re(IO"*) 


F2 


F4 


F6 


F8 


R 1 


1.307 


0.986 


0.989 


1.038 


1.001 


R2 


1.262 


0.987 


0.991 


1.017 


1.004 


R3 


0.868 


0.967 


0.997 


1.012 


0.987 


R4 


1.028 


0.978 


1.002 


1.017 


1.002 


R5 


1.177 


0.978 


1.002 


1.016 


1.003 


R6 


1.343 


0.980 


0.998 


1.016 


1.006 


K*-f 

Ka-f 


R7 


0.917 


0.998 


0.996 


1.003 


1.003 


R8 


0.935 


0.998 


0.996 


1.003 


1.003 


R9 


1.117 


0.997 


0.998 


1.002 


1.003 


RIO 


1.276 


0.997 


0.997 


1.001 


1.004 


Rll 


1.320 


0.997 


0.997 


1.007 


1.005 




0 - 270 ° 
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TABLE IX 

RESULTS OF WATER TESTS WITH SCALE 
MODEL NOZZLE TUBE 



RUN 


DATA 

POINT 


W 

lb/sec. 


NT 

in 3 / sec. 


P in 2 

lb/in 2 

gage 


A p ,. 

'tn*dis. 

Ib/in 2 


RedO* 4 ) 
d = .3l6 


^in-dis. 


WTl 


i 


0.384 


10.639 


9.175 


9.624 


2.67 


11. 187 


t 

•H 

W 

(X 

LO 

00 

r-". 

II 

B 

u 

Ph 


2 


0.495 


13.714 


14.875 


15.621 


3.44 


10.928 


3 


0.578 


16.014 


20.625 


21.642 


4.03 


11.104 


A 


0.626 


17.344 


24.200 


25.388 


4. 36 


11.105 


5 


0.656 


18.169 


26.000 


27.304 


4.56 


10.883 


6 


0.693 


19.194 


27.950 


29.411 


4.82 


10.504 


7 


0.782 


21.681 


34.150 


36.014 


5.45 


10.081 


WT2 


1 


0.254 


7.037 


3.500 


3.696 


1.77 


9.822 


P = 14.780 psi ► 

atm r 


2 


0.238 


6.594 


3.000 


3.172 


1.65 


9.600 


3 


0.382 


10.584 


8.675 


9.119 


2.66 


10. 711 


4 


0.448 


12.412 


11.950 


12.560 


3. 12 


10.728 


5 


0.521 


14.435 


16.300 


17.126 


3.63 


10.815 


6 


0.560 


15.515 


18.475 


19.430 


3.90 


10.620 


7 


0.612 


16.956 


22.000 


23.340 


4.26 


10.682 


8 


0.649 


17.981 


24.550 


25.832 


4.51 


10.513 


9 


0.690 


19.117 


27.375 


28.824 


4.80 


10.378 


10 


0.732 


20.281 


30.225 


31.856 


5.10 


10.191 


11 


0.770 


21.334 


32.950 


34.755 


5.35 


10.047 


12 


0.809 


22.414 


35. 150 


37.142 


5.64 


9.728 


13 


0.823 


22.802 


37.350 


39.412 


5.74 


9.974 


14 


0.870 


24.104 


39 . 600 


41.903 


6.06 


9.490 


15 


0.887 


24.575 


41.150 


45.545 


6.17 


9.487 
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TABLE X 

RESULTS OF WATER TESTS WITH GEOMETRICALLY SIMILAR 
MODEL OF SECTION BETWEEN PLANES W AND X 



RUN 


W 

Ib/sec 


\r 

in 3 /sec 


Re(ICf 4 ) 
d = 1.05 


f.-6 




K(.757) 
pos. B 


K(76l) 
pos. A 


KC7I4) 
pos. A 


K(.7I4) 
pos. B 


W1 


0.198 


5.488 


0.415 


1.082 


5.824 


.720 


.697 








0.313 


8.674 


0.656 


1.388 


6.886 


.768 


.727 








0.509 


14.112 


1.068 


1.167 


5.600 


.775 


.739 








0.739 


20.474 


1.549 


1.253 


5.304 


.787 


.787 








0.910 


25.203 


1.907 


1.232 


5.333 


.780 


.778 








1.007 


27.912 


2.112 


1.219 


5.375 


.774 


.771 








1.124 


31.149 


2.357 


1.219 


5.419 


.780 


.768 








1.191 


32.994 


2.496 


1.204 


5.402 


.766 


.764 








1.270 


35.203 


2.663 


1.203 


5.367 


.767 


.767 








1.372 


38.023 


2.877 


1.204 


5.576 


.755 


.753 








1.515 


41.969 


3.175 


1.208 


5.673 


.743 


.748 








1.759 


48.715 


3.686 


1.216 


5.689 


.745 


.749 








1.829 


50.682 


3.834 


1.222 


5.680 


. 746 


.751 






W2 


0.878 


24.340 


1.841 


1.224 


5.664 


. 762 


.753 








1.183 


32.769 


2.479 


1.208 


5.617 


.760 


.751 








1.344 


37.246 


2.818 


1.181 


5.718 


.752 


.736 








1.437 


39.814 


3.012 


1.178 


5.730 


.745 


.740 








1.812 


50.213 


3.799 


1.174 


5.619 


.754 


.740 








1.929 


53.435 


4.044 


1.207 


5.686 


.747 


.746 








1.946 


53.913 


4.078 


1.193 


5.768 


.742 


.737 






W3 


0.559 


15.485 


1.221 


1.487 


8.669 






.794 


.748 




0.809 


22.406 


1.767 


1.481 


8.456 






.802 


.748 




0.876 


24.285 


1.915 


1.481 


8.688 






.792 


.734 




0.988 


27.384 


2.160 


1.440 


8.418 






.793 


.746 




1.078 


29.881 


2.357 


1.462 


8.693 






.786 


.730 




1.262 


34.958 


2.757 


1.457 


8.971 






.773 


.724 




1.343 


37.222 


2.936 


1.454 


9.112 






.766 


.715 




1.440 


39.903 


3.147 


1.449 


9.035 






.768 


.718 




1.494 


41.407 


3.266 


1.452 


9.230 






.760 


.714 




1.589 


44.029 


3.473 


1.447 


9.034 






.767 


.724 




1.690 


46.828 


3.693 


1.434 


9.216 






.756 


.713 




1.754 


48.587 


3.832 


1.430 


8.698 






.778 


.733 
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TABLE XI 

OXIDIZER DISCHARGE & RETURN ASSEMBLY, INLET 
ORIFICE DATA AND EQUIVALENT FLOW RATES 

P _ - 14.743 lb/in 2 

atm 


RUN 

AVERAGE OF DATA POINTS 


R7 


R8 


R9 


R10 


Rll 


P (lb/in 2 ) 

u 


18.055 


18.943 


19.986 


21.179 


22.220 


A h Q (in H 2 0) 


0.991 


1.198 


1.448 


1.736 


1.951 


T u (°R) 


541.0 


564.5 


561.0 


554.0 


580.5 


T plen (° R > 


539.5 


559.5 


559.0 


553.5 


576.0 


Re ( 10 ' 5 ) d-4.07" 


0.917 


0.985 


1.117 


1.276 


1.320 


* 2 
w u (in ) 


.50180 


.53848 


.57723 


.61355 


. 63494 


P (lb/in 2 ) 

Z 


17.769 


18.651 


19.602 


20.750 


21.721 


P t (lb/in 2 ) 

A 


17.597 


18. 388 


19.296 


20.393 


21.330 


P t (lb/in 2 ) 

F 


15.589 


15.827 


16.114 


16.462 


16.752 


P t (lb/in 2 ) 

G 


15.464 


15.668 


15.913 


16.216 


16.463 


— * 2 
w z (in ) 


.50916 


. 54449 


. 58750 


.62596 


.64702 


w A (in Z ) 


.51414 


.55228 


.59679 


.63690 


.65888 


★ 

w p (in 2 ) 


.58037 


.64164 


.71467 


.78901 


.83895 


— * 2 
w G (in ) 


.58505 


.64817 


. 72367 


.80098 


.85364 


P t (lb/in 2 ) 

o 


18.005 


18.890 


19.894 


21.093 


22.132 



61 



TABLE XII 

LOSS COEFFICIENTS PLANES Z TO A, A TO F, 
F TO G, G TO H, H TO I, W TO X 
DISCHARGE MACH NUMBER 


RUN 


R7 


R8 


R9 


RIO 


Rll 


h-A (ln ' 4) 


.0374 


.0473 


.0451 


.0438 


.0430 


Vf (in' 4 ) 


.4316 


.4564 


.4630 


.4752 


.4944 


Vg (in " 4) 


.0237 


.0245 


.0243 


.0240 


.0244 


Vh (in ' 4) 


.0414 


.0415 


.0410 


.0407 


.0409 


P (in H O) 

I 


420.98 


425.07 


429.52 


434.54 


439.60 


Pj. (in H 2 0) 


407.00 


407.03 


406.71 


406.23 


405.83 


*1 


.22 


.25 


.28 


.31 


.34 


Vi (in " 4> 


.0500 


.0486 


.0486 


.0506 


.0493 


Vi (1 “’ 4) 


.0086 


.0071 


.0076 


.0099 


.0084 


Vx (ln ' 4> 


.4296 


.4539 


.4600 


.4702 


.4894 


P (lb/in 2 ) 

X 


15.601 


15.845 


16.138 


16.492 


16.801 


Re < 10 ' 4) d - 0.3 


2.392 


2.567 


2.913 


3.329 


3.442 


P (lb/in 2 ) 


17.338 


18.024 


18.914 


19.935 


20.752 


— * 2 

W w (in ) 


.52182 


.56342 


.60886 


.65180 


.67614 


■f =AP| w-x/ APb - c 


1.46 


1.458 


1.453 


1.447 


1.445 


S)-X <ln " 4) 


.3678 


.3908 


.3959 


.4061 


.4193 


fw-x 


9.938 


10.560 


10.698 


10.973 


11.329 
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•Combustion chamber 




NOZZLE 



FIG. 1 

ARTISTS CONCEPTION OF 
ARES PROPULSION SYSTEM. 
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SCALE MODEL OF ROCKET 
NOZZLE TUBES 
FIG. 2 
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FIG 






Orifice plates 1 Simulated pump vanes 
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IN ARES MODEL 
FIG.. 4 
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ARES MODEL 
FIG. 7 
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VIEW OF RETURN ANNULUS 
FIG. 10 






VIEW GF SUPPLY PIPE 
A: ID PLENUM 2 
FI I 






A 




A - Pressure upstream 
B - Pressure downstream 
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Fig 12 Schematic of sharp edge Orifice 
and flow nozzle. 
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SCHEMATIC OF FLOW MEASURING 
INSTALLATION 
FIG. 13 
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USED AT MODEL INLET 
' FIG. 15 




SCHEMATIC OF INLET TOTAL PRESSURE 
PROBE INSTALLATION 
FIG. 16 
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IN EITHER RADIAL OR 
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FLOW 




SCHEMATIC OF INJECTOR TUBE AND TOTAL 
PRESSURE PROBE ALIGNMENT 
FIG. 18 
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VIEW OF A FLOW CHANNEL 
FIG. 22 
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INLET TO DISCHARGE LOSS COEFFICIENT OF FLOW DEVICE 

VS REYNOLDS NO. 

FIG. 24 




OF FLOW DEVICE 
FIG. 25 
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FIG. 30 
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PLANE D TO PLANE X. 
FIG. 31 
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FIG. 32 

FRICTION FACTORS FOR FLOW IN PIPES 
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FIG. 33 

PERIPHERAL DISTRIBUTION OF LOSS COEFFICIENT 
PLANE G, RUNS Rl THRU R6. 




VIEW OF INTERSECTION OF PASSAGES FROM 
PLANE G TO RETURN ANNULUS 
FIG. 34 




Re x ICT® 

AVERAGE LOSS COEFFICIENTS RETURN SUB-ASSEMBLY 

VS REYNOLDS NO. 

FIG. 35 
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FIG. 37 
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DISCHARGE VELOCITY DISTRIBUTION RUN R4 

FIG. 38 
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DISCHARGE VELOCITY DISTRIBUTION RUN R5 

FIG. 39 
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DISCHARGE VELOCITY DISTRIBUTION RUN R6 

FIG. 40 



FIG. 41 



ISIS 
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FIG. 42 FLOW IN RETURN ANNULUS RUN R4 
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FIG. 42 (con't) FLOW IN RETURN ANNULUS RUN R4 
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FIG. 42 (con't) FLOW IN RETURN ANNULUS RUN R4 
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FIG. 42 (con't) FLOW IN RETURN ANNULUS RUN R4 
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FIG. 43 FLOW IN RETURN ANNULUS RUN R5 
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FIG. A3 (con't) FLOW IN RETURN ANNULUS RUN R5 




FIG. 43 (con't) FLOW IN RETURN ANNULUS RUN R5 
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FIG. 43 (con't) FLOW IN RETURN ANNULUS RUN R5 
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NOZZLE TUBE LOSS COEFFICIENT VS REYNOLDS NO. 
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FRICTION FACTOR VS REYNOLDS NO. FOR RUNS Wl AND W2 

FIG. 45 
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FRICTION FACTOR VS REYNOLDS NO. FOR RUN W3 

FIG. 46 
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PERIPHERAL DISTRIBUTION OF LOSS COEFFICIENT 
PLANE G, RUNS R7 THRU Rll. 
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APPENDIX A 



DERIVATION OF FLOW MEASURING EQUATION 



The pressure drop P^ - P^ across the flow nozzle is referred to 

the absolute pressure P upstream and is to enter calculations in the 

Pu-Pd 

dimensionless form — — — . Thus any consistent unit of pressure may be 



Pu 



used. 



The equation for flow rate w (lb^/sec) is 



w 



- *■ X Aj f t f r Hjz f ^ (P u - P d ) 



(Al) 



where: 

oL - nozzle discharge coefficient (independent of Re) 

- expansion factor to account for compressibility effects 

2 2 

A^ - area based on nozzle diameter (in. ) 3 7. 17263 in. 
f fc - area multiplier to account for the thermal expansion of the 
nozzle 

f^ - Reynolds number correction factor 

g - 32.17 ft/sec^ 

Pu - upstream density (lb /ft ) 

m 

2 

Pu - pressure upstream (lb/in. ) 

2 

Pd - pressure downstream (lb/in. ) 

Using the equation of state for a perfect gas 



P - p * Cm-*') 

V u 

* 



(A2) 



1943, 



^VDI Durchfluss 
p. 5. 



Me ssre geln. 



DIN 1952, VDI Verlag, Berlin, 
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where R = 53.35 ft-lb/lb^R, and is the upstream temperature (°R) , 
equation (Al) can be rearranged to give the equivalent flow rate 



vs/ 



^ \[E = Y, A d f t f 



U 






~ (if?) (A3) 



The dimensionless coefficients must be investigated for the noz- 

zle with diameter ratio ^ * 0.7434 and m = ^ = 0.552 where = d/D. 

2 

The discharge coefficient oC for m » 0.552 is found to be 1.114. 

3 

The density coefficient is obtained from the ASME standards. 

For air with ^ = c /c = 1.4 and ■ .743 using x = ~r~ there is 

p v V Pu 



Y. - 



2.433 



0- X. 



/.+za i>.Z86 

) /~ ] 
( / - 0.305 ( I - K) / " f2a J 



(A4) 



L. 

For stainless steel the area multiplier f^ can be expressed by 



■f = /, 0 + O.OOtS f 

i I I 



t>. - 60 



o o 



(A5) 



The design of the nozzle is such that f^ = 1 for all pressure ratios. 
Combining the above quantities equation (A3) becomes 




II. 3 




(A6) 



2 

Ibid . Appendix 6. 

3 

Supplement to ASME Power Test Codes, Chap. 4, Part 5, ASME, 
New York, 1959, p. 76. 

4 

Stearnes, R. F. et. al. Flow Measurements with Orifice Meters , 
Van Nostrand Co., Inc., New York, 1951, p. 257. 
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APPENDIX B 



DATA REDUCTION FOR AIR AND WATER TESTS WITH FLOW CHANNEL 



B. 1 Air Tests 

A total of seven data points was observed. The quantities mea- 
sured were: 

H - pressure upstream of flow measuring orifice against atmo- 
sphere (in-Hg) 

Ah Q - change in static pressure across orifice (in^H^O) 

t^ - temperature upstream of orifice (°F) 

t plen ~ temperature in plenum 2 (°F) 

H - manometer reference (cm H o 0) 
atm l 

H^ n - inlet total pressure at Kiel probe against atmosphere 
(cm H 2 0) 

H^ s - discharge total pressure against atmosphere (cm ^0) 

h 30.01 in. Hg 

atm 

The equivalent flow rate w^ is given by equation (A6). It is 
convenient to use the model inlet as a reference condition. Then the 
equivalent flow rate is 



w 



w 



i f\ 



IA 



t. 



It 

T 



_LH 






(Bl) 



IA 



where P y = P,. = 0.49112 j H y + 30. Oil (lb/in 2 ) 

u * * 

P, = P, = 0.49112 ( 0.02894 (H - H 4 ) + 30.01 \ (lb/in 2 ) 
in t. I v atm in' ) x 



Reynolds number Re based on inlet diameter d is 



ie - 



Vg) ^ 

/< * 



4 w 



IT 



* % 



(B2) 
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Equation (B2) can be rewritten as 



* 



Re -- 



w 

in 




4 ( I 2 ^ 
Wd / i / ** 



(B3) 



Air temperature for these tests varied between 85 and 95°F. Introducing 

-7 2 

a representative viscosity of 3.9 x 10 lb-sec/ft , in addition to 
other non-varying quantities, into equation (B3) gives 



Re * w i* (4.7£e, * /o s ) 




The loss coefficient is determined by 



(B4) 



( in 4 ) < B5 > 

where P t = 0.49112 j 0.02894 (H - H ) + 30.0l\ (lb/in 2 ) 

tj, ( atmt,. j 

dis v dis 

B . 2 Water Tests 

A total of four runs was made, each consisting of several data 
points. The runs were numbered RMW 1, etc. The quantities measured 
were: 



in-cUs 







/A 



2 

g - inlet static pressure (lb/in gage) 
W - weight of water (lb) 
t - elapsed time (sec) 

Volume flow rate V** is obtained by 



V - 



Vv ( 100 ) 

3.6093 



( in* / S«c) 



(B6) 



Re based on inlet diameter is 

Re = =. 'f. 

M y ^4 d V 



(B7) 
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where V <* kinematic viscosity, for water at 72°F * 152.257 in^/sec 



Equation (B7) is reduced to 
Rd - \P ^ 2.39. 76 4-5 

Because the inlet and discharge areas are equal 



(B8) 



- AP. - P. 

t. m-dis in g 

in-dis 



The loss coefficient is determined by 



k- 



i(\- di* 



P Y. 

'L 

\v 2 



/ 3 .9 4-O-b P. 



M 



w 



(irt"*) (B9) 



where Y = specific weight of water at 72°F. 
Li 
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APPENDIX C 



DATA REDUCTION FOR TESTS WITH ARES MODEL 



C. 1 General 

Methods of data reduction used to obtain loss coefficients and 
other meaningful quantities from tests on the ARES model (as described 
in Section 5 of the main body), are discussed in the following Sections. 
They are listed in the same sequence as used for the performance of the 
tests . 

C.2 Oxidizer Return Sub-assembly 

The quantities necessary to obtain mass flow rate are identical to 

* 

those described in Appendix B. The equivalent flow rate w^ is given by 
equation (A6) . In addition the following quantities were recorded 
Hpien ~ pressure in plenum 2 (cm ^0) 

H - - atmospheric reference, manometer bank 1 (cm H o 0) 
a tm 1 Z 

H, - pressure at Kiel probe (cm H o 0) 
in Z 

- pressure in plane D (cm 1^0) 

- static pressure at selected locations in annulus (in. 1^0) 

H an - total pressure at selected locations in annulus (in.H^O) 

H _ 0 - atmospheric reference, manometer bank 3 (in.H o 0) 

a tm 3 2 

hj - static pressure in plane I (in. ^0) 

Hj - total pressure in plane I (in.I^O) 
tpien “ temperature in plenum 2 (°F) 

From the polaroid pictures the following quantities were read 

hg - static pressure in plane E at holes 7, 13, 20, 26, 33, 39, 

46, 52 (cm 1^0) 

h - static pressure in plane F at holes 1, 9, 17, 25 (cm H o 0) 
r Z 
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H, - total pressure in plane G between holes 32-1, 4-5, 9-10, 
G 



13-14, 17-18, 20-21, 25-26, 29-3 (cm H 2 0) 

h - static pressure in plane H above discharge holes 11, 14, 18, 
H 

22, 26, 30, 31, 34, 38, 42, 46, 50 (cm H 2 0) 



Re * w u* Jk. j 

/nr it m . 



w* P * ‘ > - 0,ofc * (Cl) 

“ fZ A 



The loss coefficient k for flow between each plane is referred to 
an inlet condition. The equivalent flow rate must then be referred to 
conditions at intermediate planes. For example^an average equivalent 
flow rate at plane F is 



W 



w\ 



IA 



Tr 



(C2) 



P, V T 

where ? t - 0.49112 ^ 0.02894 (H atn x + P atm (lb/in 2 ) 

F 1 

T^ - temperature in plane D assumed to be the same as that in 
plenum 2 



P and T - determined as in equation (Bl). 
u u 

___ ^ n ^ 

In a similar fashion the equivalent flow rates w^ , w^ are determined. 



A geometrically similar model of a tube section from plane D to 
plane X was tested by Vavra. Data from these tests plotted in Figs. 29 
and 30 include the orifice discharge coefficient K, and the loss in 
total pressure ^ from plane D to plane X divided by the static 

pressure change across the orifice. For each orifice preceding a mea- 
suring tap in plane E, an equivalent flow rate is found by 
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* 



w 



o 



« Y, A„ 





(C3) 



2 

where A, =0. 0363 in (area of orifice opening). 



P and P are determined 
E D 

® * c /c *1.4 from 
P v 



a 1.0 - 0.36 

i 



Equation (C3) becomes 



* 

W - 0,0$ Z0+ l< Y, 

0. 

U 



as was P above. For air with 
F 




(C4) 



(C5) 



Equation (B4) is modified to find Re referred to tube inlet diameter 



Re « 



P / 7 \ 

o ( <3.2>Z 3 x io ) 



(C6) 



This relation contains a representative value of viscosity for the tem- 
perature range encountered in the tests. An iteration is carried out 
using equations (C5) and (C6). A value of K is assumed and equation 



(C5) is solved. Using this value of w , Re is calculated from equa- 

D i 



tion (C6) . A corrected value of K is then read from Fig. 29 and the 
process continued. Generally, one iteration is required only. 

The peripheral flow distribution at the orifice plate is found 



from 



F. 



(£z) 



* 



w 



o » 

4 




— # 



w 



0 



(C7) 
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Tht total pressure P is found from 

C X 



=> = P 

-t_ ° 



- { ( ' P ) 
Jd-x V D 1 ' 



(C8) 



where £ is read from Fig. 30 for the Re found with equation (C6). It 

is possible also to estimate P fc by using equation (C7) with £ = 0.48, 

E 

a standard value for sharp edge orifices with ^ = 0.72. 

The flow area in plane E is 3.676 sq. in. It is possible to find 

an average area restriction factor . The equation of continuity for 

E 

average conditions is 



v 



w 



e a <£ 



(C9) 



Also 



_ 2, 

f \ 



- j \ h - p \ 



(CIO) 



Combining equations (C9) and (CIO) gives 



W 






p-p| 



(Cll) 



Using the equation of state for a perfect gas P 
becomes 



RT, equation (Cll) 



w trr 

P 2 q. 






tv 



(C12) 
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Referring equation (C12) to static conditions at plane E and solving 
for 

^2 






, 1*1) 






(Cl 3) 



where : 



V, 



w 






An average Mach number M is determined from the pressure ratio 

P /P . Assuming no change in total temperature across the orifice, 

E fc E 

i.e.; T 3 T , it is possible to find the temperature ratio needed to 
E D 

solve equation (C13). 

With the equivalent flow rates and equation (12), average loss 

coefficients lc , k , k , k and k are determined. For example 
D - F F~G G~ri G~ 1 D ~X ' 

the average loss coefficient from plane F to plane G is 






P-G 



G_ / 2 

(Z r *) 



.(.V 4 ) 



(C14) 



At a particular location in plane I the discharge velocity is pro- 
portional to the square root of the dynamic pressure. 



V.*. [h 

1 L 






ff\ 3 



H c'] ‘ [ (c i5) 



At each location the velocity relative to the average is then 

V; / _ _ 



vy 77 ZfF* 



(C16) 
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C. 3 Scale Model Nozzle Tubes 

The following quantities were recorded 
t* - elapsed time (sec) 

W t - weight of water collected (lb) 

2 

P, - inlet static pressure (Ib/in gage) 
in g 

t - water temperature (°F) 

2 

A, - inlet area (0.07843 in ) 
in 

2 

A^ g - discharge area (0.11341 in ) 
h^ ” barometric pressure (in Hg) 

Flow rate w was determined by dividing the weight of water by the 
elapsed time. More useful for calculations is the volume flow rate 



\F = 



w 



(lOO ) 



( in* / sec.) 



3. 6 <3?2> 



(C17) 



Velocities at inlet and discharge can be written as 

V = / A (ci8) 

Assuming atmospheric pressure at the discharge the loss in total 
pressure through the tube is 



P ‘d, 5 * I Pi "~ * 2 (V ' K ^ 



(C19) 



In all tests the water temperature varied only from 55° to 59°F. 



For 



t ■ 57 F equation (C19) becomes 



a n 



- 3 



^in-c/i's 



P;, , y, +■ ( x /o' ) \f 2 



(C20) 



Re referred to an inlet diameter of 0.316 inches is obtained from 

-5 . 2 



equation (B7) with V - 160.279 x 10 in /sec 






\f * /o 

sTa 



(C21) 
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Using the inlet as reference, the loss coefficient 




* 

* /n - di 



* i n - d i $ 



(z 







z 



(C22) 



where 




-4 

0. ?3S x /O 



z 

Jb~ se c 



/V) 



4 



C.4 Geometrically Similar Model of Section from Plane W to Plane X 

Figure 28 shows 6 manometer tubes in a schematic manner. Actually, 
Stations 2 through 6 have 2 taps each so that h^ (in-H^O) is an average 
pressure at Station 2. In addition to pressures, the mass flow rate 
and water temperature were recorded as described in sub-section C.3. 
Volume flow rate is determined as in the former tests. Reynolds number 
referred to an inlet diameter of 1.05 in. is 



\F * < o 4 

13 .2 IS 



(C23) 



for average tempera tue conditions. 

It is assumed that taps 3 and 5 are placed such that the flow has 
the same area as the orifice opening. The error in this assumption is 
absorbed in the flow coefficient K. 

Across an orifice the equation of continuity gives 

A d V dl = (C24) 

where d - station ahead of orifice 

or - station behind the orifice 

Expressing terms in in. of liquid head, Bernouilli's equation, ignoring 
losses, gives 
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z 



(C25) 






k 



Combining equations (C24) and (C25) 





( Kj- K^~| 

'-P* 



(C26) 



If losses and other factors are included in a coefficient K, the volume 
flow rate can be written as 



v - y e = k v- or 

Equations (C26) and (C27) are combined to give 



(C27) 



K = 








(C28) 



Three orifice sizes were tested. Table C-I lists factors for each 



orifice . 



TABLE C-I 



GEOMETRIC 


FACTORS 


FOR VARIOUS 


ORIFICES 


DIAMETER 


6 


1-fc 4 


A 




V 


T 


or 


0.757 


0.7209 


0.7298 


0.45007 


0.761 


0.7248 


0.7241 


0.45484 


0.714 


0.6800 


0.7862 


0.40039 



Using Station 6 as an example and assuming stagnation conditions 



at Station 1, 



the head loss 



«L 



is written as 



‘ ^4 = H , ~ \ ' V JS- 






/ * 2 



1 



(C29) 



Compared with otfier losses, the loss through the inlet orifice is 
negligible. A loss coefficient from the beginning of the 1.05 in. pipe 
to Station 6 is then 




(m 







Introducing the volume flow rate 



(C30) 



- * 78 - 8 ’ 6 

\F Z 



(C31) 



C. 5 Discharge and Return Assembly 

In addition to the quantities measured in the return sub-assembly 

the following pressures were read 

h Q - static pressure at entrance of model (cm 1^0) 

H q - total pressure at entrance of model (cm 1^0) 

h 0 - atmospheric reference, manometer bank 2 (cm H„0) 
atm 2 2 

h z - static pressure in plane A, 6 equally spaced taps (cm ^0) 

- total pressure in plane Z (cm 1^0) 
h^ - static pressure in plane A at holes 1, 14, 27, 40 (cm 1^0) 

- total pressure in plane A between holes 7-8, 20-21, 33-34, 
46-47 (cm H^O) 



h 



B 



h - 
c 



H plen 



static pressure in plane B at holes 7, 13, 20, 26, 33, 39, 
46, 52 (cm H^O) 

static pressure in plane C at holes 7, 13, 20, 26, 33, 39, 
46, 52 (cm H^O) 

- pressure in plenum 2 referenced to atmosphere (cm Hg) 
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Equivalent flow rates found as in Section C.2 include w , w , 

n U 2 

-- ^ -- «n< 

w , w_ , and w . These flow rates are used to calculate the loss 
A F G 

coefficients ^_ A> k A _ p , k^, k^, and k^. 



Using Re with the 4 in. diameter of the supply pipe as character- 
istic length, the argument, k^_^ is read from Fig. 35. Total pressure 
at plane X is then given by 



X 




(C32) 



where : 



k 






=. Ir - L 

A-F K X-F 



An average Re of the flow from plane W to plane X can be written as 



where: 



Re - ^ 



— * ft, (^.02364-) 



"V 



V 



(C33) 







(C34) 



and 



p t 



\r 



- { (* p *..) ♦ p* 



(C35) 



An iteration of equations (C33), (C34), and (C35) is carried out by 
assuming an initial value of ^ from Fig. 46. The process converges 
rapidly because has a limited range of values. 

With the computed values of P^_ and w T *. the coefficient „ is 

S? w 

calculated. From equation (12) 

^ W -X “ A ( C36 ) 



W - A 

Z 4 l 
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